Exploratory spectral analysis of hydrological times series     

A. Ian McLeod  K. W. Hipel 《Stochastic Environmental Research and Risk Assessment (SERRA)》1995,9(3):171-205

Current methods of estimation of the univariate spectral density are reviewed and some improvements are made. It is suggested that spectral analysis may perhaps be best thought of as another exploratory data analysis (EDA) tool which complements, rather than competes with, the popular ARMA model building approach. A new diagnostic check for ARMA model adequacy based on the nonparametric spectral density is introduced. Additionally, two new algorithms for fast computation of the autoregressive spectral density function are presented. For improving interpretation of results, a new style of plotting the spectral density function is suggested. Exploratory spectral analyses of a number of hydrological time series are performed and some interesting periodicities are suggested for further investigation. The application of spectral analysis to determine the possible existence of long memory in natural time series is discussed with respect to long riverflow, treering and mud varve series. Moreover, a comparison of the estimated spectral densities suggests the ARMA models fitted previously to these datasets adequately describe the low frequency component. Finally, the software and data used in this paper are available by anonymous ftp from fisher.stats.uwo.ca.  相似文献   

基于误差准则和循环迭代的时移地震匹配滤波方法   总被引：13，自引：3，他引：13       下载免费PDF全文

金龙  陈小宏  李景叶 《地球物理学报》2005,48(3):698-703

分析了用于时移地震互均衡处理中的传统匹配滤波算法适用范围，结果表明该方法在滤波算子的反算子为最小相位时效果最好.对该方法难以解决的混合相位问题，推导出通用公式，并提出基于最小平方误差准则和循环迭代的求解方法.理论和实际数据都验证了该方法较已有方法的优越性.  相似文献   

Transient deformation following the 30 January 1997 dike intrusion at Kīlauea volcano,Hawai’i     

Emily K. Desmarais  Paul Segall 《Bulletin of Volcanology》2007,69(4):353-363

On 30 January 1997 an intrusion on Kīlauea volcano opened a new fissure within the East Rift Zone (ERZ) at Nāpau Crater, 3 km uprift from the ongoing eruptions at Pu’u ’Ō’ō. The fissure eruption lasted 22 h and opened a 5.1 km long, nearly vertical dike 1.9 m, extending from the surface to a depth of 2.4 km (Owen et al. 2000b). During the eruption, the lava pond at Pu’u ’Ō’ō drained, and eruptions ceased there. Pu’u ’Ō’ō eventually refilled in late February and eruptions resumed there on 28 March 1997. Continuous GPS data show a large transient following the 30 January 1997 dike intrusion. After lengthening 40 cm during the initial eruption, the baseline between two stations spanning the ERZ lengthened an additional 10 cm over the following 6 months. A coastal station KAEP also exhibited transient deformation, as it continued to move southward (5 cm) over the same 6-month period. The baseline between two stations spanning Kīlauea’s summit caldera contracted sharply during the eruption, but gradually recovered to slightly longer than its previous length 2 months after the intrusion. We use the extended network inversion filter (McGuire and Segall 2003) to invert continuous GPS data for volume change of a spherical pressure source under Kīlauea’s summit, opening distribution on a nearly vertical dike in the ERZ and potential slip on a decollement 9 km beneath the south flank. Following the 30 January intrusion, rift extension continued below the initial dike intrusion for the duration of the transient. Decollement slip, regardless of its assumed depth, is not required to fit the data. The modeled transient summit reinflation and rift opening patterns under Nāpau crater coincide with changes in observed behavior of Pu’u ’Ō’ō’s lava pond. Rift opening accelerated while Pu’u ’Ō’ō eruptions paused and began to decelerate after the lava pond reappeared nearly a month after the Nāpau eruption. The transient deformation is interpreted as resulting from shallow accommodation of the new dike volume.  相似文献   

Walsh列率域中多维分形模型与GIS环境下地球物理信号处理   总被引：4，自引：0，他引：4  

李庆谋  成秋明  刘少华 《地球物理学报》2007,50(6):1884-1893

地球物理信号通常在多个尺度段表现尺度不变性，这些不变性起因于不同的地质、地球物理或成矿过程的自相似性. 利用这种在多个尺度段的尺度不变性可以设计多维分形滤波器，滤波所得信号表征了具尺度不变性的地质地球物理或成矿过程，可以用于成矿预测或环境评价. 本文研究了Walsh变换列率空间地球物理信号的列率功率谱密度与列率之间的分形与多维分形关系, 试验证实了大洋钻探、石油以及煤系地层地球物理测井资料在Walsh域的多维分形性质，提出了用于分解地球物理场,提取有用信号并用于矿产资源勘探或环境评价的多维分形W-A模型. 利用波列率域中的多维分形关系构造了W A图解(W-A Plot). 借助W-A图解可以确定最小平方误差（LS）意义下Walsh功率谱变化的不同自相似性的频率分界点,从而用于设计W-A分形滤波器. 这种滤波器可将地球物理场分解成具有不同自相似性的局部场,并且保留原场的各向异性结构. 与通常使用的基于Fourier变换的滤波技术相比,W-A模型具有许多优点：W-A适用于检测地球物理信号中的突变、线性、环状、局部与纹理结构等弱信号. 同时，由于Walsh变换中只有简单的变号(加法与减法)，其计算速度远快于建立在复数乘法之上的Fourier变换，所以W-A计算速度远快于Fourier域的滤波方法，可以用于地球物理信号的现场实时处理. 用加拿大Nova-cotia省西南地区的布格重力异常进行了W-A方法的试算，处理结果反映了地质、矿产分布规律，能够很好地进行矿产预测.  相似文献   

对重力势二阶垂直偏导数的一种新的取样方法     

AHMED  A.  KHAMIES  许军才 《火山地质与矿产》2007,28(4):242-248

从目前的研究来看，二阶垂直导数（X和Y两个方向）的求解方法较简单，它采用的五个因子的数字滤波的最小二乘法来优化数据。但是在较小区域面积时（对工程和环境应用中的重力数据），重力异常很容易放大和被区域重力场掩盖，我们认为产生的缘由是表示重力场的x和y方向上的二次多项式引起的。因此，为了逼近重力水平方向的偏导可以设计一个好的滤波器来拟合在xy平面的重力场值，再得出想要的重力垂直方向偏导，文中最后对提出的方法在前人经典模型上和野外实践进行了验证。  相似文献   

基于启发式图搜索的遥感影像道路半自动提取   总被引：10，自引：2，他引：10  

李燕  余旭初 《测绘学院学报》2002,19(4):268-271

启发式图搜索法用于线状目标识别的原理是：用图结构表示边缘点和边缘段，根据启发函数计算顶点权值，在图的路径上建立相应的代价函数，通过在图中搜索对应的最小代价的通道以找到最优路径。图搜索法是一种全局最优方法，它在受噪声影响较大时效果仍然较好。文中使用了启发式图搜索法（A*算法）实现了道路的半自动跟踪。它的基本思路是：首先利用自适应平滑滤波算子进行道路信息增强，然后对传统的道路数学模型进行了进一步的扩展，突出了对道路几何特性和辐射特性的描述，并依此构建图搜索的代价函数，实现了基于启发式图搜索法A*算法的道路半自动跟踪。经实验证明，该方法进行遥感影像的道路半自动提取效果较好。  相似文献   

华北平原深层淡水在开采条件下接受上覆咸水越流补给——以天津平原为例   总被引：1，自引：5，他引：1  

王家兵 《水文地质工程地质》2002,29(6):35-37

分析物探测井曲线和地下水位动态发现，天津平原深层孔隙地下淡水在开采条件下接受上覆咸水的越流补给，而且咸水底界未大量向下位移。本文讨论了咸质含水组底部钙化粘土微观粘粒与水溶液的相互作用，提出了“滤盐层”的概念。  相似文献   

多级汇水盆地的自动提取   总被引：1，自引：0，他引：1  

余桂贤  孙光耀  郑凯梅  吴磊  余生晨 《物探化探计算技术》2002,24(4):349-351

利用多级汇水盆地在频率城（或波数域）中的频率差异，通过滤波方法逐级分离，并提取它们，实验证明，该方法效果较好。  相似文献   

SINS／GPS级联组合滤波新算法     

刘广军  吴晓平 《测绘学院学报》2001,18(3):176-178

针对级联式SINS／GPS导航系统中组合滤波器的测量噪声与系统噪声的相关问题，提出了一种相关程度未知条件下的估计量最优融合算法，并由此得到了一个顾及噪声相关性的的卡尔曼滤波新算法；通过对一套SINS／GPS数据的计算，证明了新算法是成功的。  相似文献   

数字滤波器的频谱窗和噪声压缩度   总被引：2，自引：0，他引：2  

彭军还 《测绘学院学报》2001,18(4):246-249

根据傅里叶理论导出了数字滤波器的频谱窗，根据随机过程理论导出了数字滤波的噪声压缩度。通过对二项式系数数字滤波器和多项式系数数字滤波器和频谱窗和噪声压缩度的分析，发现前者对噪声具有较强的压缩作用，同时能保持信号的凹凸性，是一种优于后者的数字滤波器。  相似文献